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Abstract 

The search strategy or the discovery of new effects for heavy neutrinos often rely on their different 
decay channels to detectable particles. In particular in this work we study the decay of a Majorana 
neutrino with interactions obtained from an effective general theory modeling new physics at the 
scale A. The results obtained are general because they are based in an effective theory and not 
in specific models. We are interested in relatively light heavy Majorana neutrinos, with masses 
lower than the W mass {rriM < mw)- This mass range simplifies the study by reducing the 
possible decay modes. Moreover, we found that for A ~ 1 TeV, the neutrino plus photon channel 
could account for different observations: we analyze the potentiality of the studied interactions to 
explain some neutrino-related problems like the MiniBooNE and SHALON anomalies. We show 
in different figures the dominant branching ratios and the decay length of the Majorana neutrino 
in this approach. This kind of heavy neutral leptons could be searched for in the LHC with the 
use of displaced vertices techniques. 
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I. INTRODUCTION 


One of the most spectacular new results in high energy physics is the discovery of neutrino 
oscillations, indicating that they are not massless. The neutrinos can be of two different 
types: either Dirac or Majorana particles. Dirac fermions have distinct particle and antipar¬ 
ticle degrees of freedom while Majorana fermions make no such distinction and have half 
as many degrees of freedom [1]. In this conditions fermions with conserved charges (color, 
electric charge, lepton number,...) must be of Dirac type, while fermions without conserved 
charges may be of either type. New still undetected neutrinos could have large masses and 
be of either type. If heavy neutrinos (N) do exist, present and future experiments would 
offer the possibility of establishing their nature. The production of Majorana neutrinos via 
e“ 7 , 77 , e~P and hadronic collisions have been extensively investigated |2HI9]. 

A very known scenario for the study of Majorana neutrinos is the seesaw mechanism 
[20] , requiring the existence of at least one type of heavy right-handed Majorana neutrino. 
As indicated in (TH], the parameters determining the interaction of the heavy Majorana 
neutrino N with the standard particles turn out to be very small, indicating the need for a 
new approach involving physics beyond the typical seesaw scenarios. 

In this work we study the decay modes of a relatively light heavy Majorana neutrino in 
the context of a general effective framework. We focus in a mass interval below the standard 
massive vector bosons mass (m^r < m\y) as this reduces the possible decay channels, letting 
us concentrate on the phenomenology of the neutrino plus photon mode. This heavy neutrino 
decay channel has been introduced as a possible answer to some experimental puzzles, like 
the MiniBooNE l22] and SHALON |23] anomalies, considering sterile heavy neutrinos 
created by neutral current interactions and decaying radiatively due a transition magnetic 
moment [23]. We revisit here the mentioned anomalies in the light of an effective Lagrangian 
description for the heavy Majorana neutrino decays. 


The paper is organized as follows: in subsections IA, IB and IC we describe the model- 
independent effective approach, and show our analytic and numerical results for the Majo¬ 
rana neutrino decay widths and branching ratios, and discuss the existing bounds on the 
effective couplings. In section [IT] we explore the potentiality of the effective approach to 
explain the MiniBooNE and SHALON anomalies, showing our results for the Majorana 


neutrino lifetime and decay length. We present our conclusions in section III 
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A. Effective Lagrangian 


As it was explained in |Tn], the presence of Majorana nentrinos would be a signal of 
physics beyond the minimal seesaw mechanism, and thus their interactions would be best 
described in a model-independent effective approach. It is possible to parametrize the effects 






FIG. 1: Decay modes for Majorana neutrinos with mv < mw 

of new physics beyond the standard model by a set of effective operators O constructed with 
the standard model and the Majorana neutrino helds and satisfying the Standard Model 
SU{2)l®U{1)y gauge symmetry [19]. The effect of these operators is suppressed by inverse 
powers of the new physics scale A, for which we take the value A = 1 TeV. The total 
Lagrangian is organized as follows: 

OO . 

( 1 ) 

n=6 i 

For the considered operators we follow [19] starting with a rather general effective Lagrangian 
density for the interaction of right handed Majorana neutrinos N with bosons, leptons and 
quarks. The following are dimension 6 operators and could be generated at tree level in 
the unknown fundamental ultraviolet theory. The hrst subset includes operators with scalar 
and vector bosons (SVB), 

Oln^ = = i{f eD^ct>){Nrk) ( 2 ) 
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and a second snbset inclndes the baryon-nnmber conserving 4-ferniion contact terms: 


OduNe = , Ojnn = (3) 

OiNLe = {LiN)e{Lili) , OiNQd = {LiN)e{Qidi), 

OquNL = {QiUi){N Li) , OqNLd = {QiN)^{Lidi), 

Oln = \NL,\^ 

where k, Ui, di and Li, Qi denote, for the family labeled i, the right handed SU{2) singlet 
and the left-handed SU (2) donblets, respectively. In addition, there are operators generated 
at one-loop level in the nnderlying fnll theory whose coefficients are natnrally snppressed by 
a factor l/lbvr^ [T^ 125] : 

Onb = {La^'^N)4>B^,, Onw = N)^Wl,, 

Odn = {LD^N)Dy, Odn = {D^LN)Dy . (4) 

Taking the scalar donblet after spontaneons symmetry breaking as 0 = ^, with h being 

the Higgs held, the operators listed in ([^ contribnte to the effective Lagrangian 

laj’ f i^L.iNR H.iNR ^ hhh) 

-az (^-(NrYNr) + 

+ (NrYNr) j ) 

- + /i.cj . (5) 

The 4-fermion Lagrangian can be written ([^: 

= — ^a^Y^dji^i'j^UR^iNji'j^lji^i + a^yjR^i'y^lR^iNR'y ij,NR + a^y^Li'y^ LiNR'y r+ 
OivluR^i'y^UR^iNR'jfj^NR + ay^dR^i'j^dR^iNR'y^NR + ay.Qi'^^QiNR'-^^NR+ 
<^sli^L,iNReL,ilR,i — eL,iNRUL,dR,i) + 0'gli'dL,iUR^iNiyL,i + dL,iUR^iNeL,i) + 
(^sli^L,iNRdL,idR^i — CL^iNRUi^idR^i) - 1 - agliuR^iNRei^idR^i — dL,iNRyL,idR^i) + 

iL,iNR + NReL,ieL,iNR) + h.c.'^ (6) 
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In Eqs. 0 and ^ a sum over the family index i is understood, and the constants a^Q are 
associated to specific operators: 


az — 
(i) 

= 

h) 

= 


OiNNtj), — OL 


» _ ^,h) 

LN4>^ 


(i) _ (h 

“w — ^ATec/.’ 


(d _ (d 

*^Vo “ ^duNe’ 


(d _ (d 

“ ^eNNi 


„(d 


„(d 


„(d 


.,(d 


„(d 


., _ ., .. _ ., (d _ (d ., _ . 

^LNNy ~ ^uNNy *^^4 ~ ^dNN^ ~ ^QNNi *^So ~ ^LNe^ 


„(d 


„(d 


„(d 


„(d 


„(d 


„(d 


„(d 


„(d 


„(d 


^QuNLi ^S 2 ~ ^LNQdi ^83 ~ ^QNLdi ^ 84 , ~ ^LN' 


(7) 


For the case of the one-loop generated operators in Q, we have the effective Lagrangian: 


(-iV2vcwPl^'^ ^L,^a^''NR A, + iV2vswPi^'’ ^L,^a^''NR Z,+ 
-iV2cwPl,^^ H,^a^'^NR A,h + z Nr Z,h) 

+mwPl4^^ Tl,Nr W-^ + 1^,Nr W-^h + Vl,Nr h 

[i^PlvcwPf^ vl^cj^'^Nr Z, + iV2vswP!4^^ dl, 

+t2V2mw y^^(T^''NR IF+IF- +zy2nP^'^) Tl,(t^'^Nr IF" 

+iAmwCw lL,iO'^‘'NR Zy + iAmwSw Il^o'^’^Nr Ay 

+zV2Pl^^^ ^L,^a^^NR W;h + i2gcw lL,^a^''NR WyZ^h 
+i2gsw h,iO'^‘'NR W~A^h + i\/2cwPjf^ Z^h 

+iV2swPl^'^ VL,iCT>^''NR A^h + iV2g PL,ia^'^NR W+W-h}j 


a 


,(d 


A2 

1 


(^^7’ ~ N'’) Z^h 

^IpC'.jcpW h - 228A w-^w; - ^ n,.iV„ Z^Z" 

v2 V y/2v 

ly^iNR z,z^h - ur^^Nr w+w-^h 

+mwPpW->^ Il,Nr + - Pf ))IF-'^ In.iVK h 

+emw lL,iNR W ^A^ + emzsw Il^Nr W ^Z^ 

- z^W-^h + Il^^Nr A,W-^h\ + h.c. 

V2n / 

where P*^“^ is the 4-moment of the incoming a-particle and a sum over the family 
i is understood again. The constants 


( 8 ) 

index 


af. with j = I, 


O Q 
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operators: 






„(0 _ Ji) 


„(*) _ ^,(*) 


*^Li ~ *^7V_B! '^L2 ~ '^DAT! ~ ^NW! ^La ~ ^DN' 


(9) 


B. Decay Widths 


We have calculated the decay channels for a Majorana neutrino with mass lower than 
the standard model vector bosons {rriN < mw)- This range allows for the decay to fermions 
(excepting the top quark) and to photons. The contributing decay modes are schematically 
shown in Figj^ 

We now present the partial decay widths of a heavy Majorana neutrino N decaying to 
three fermions. They were calculated using the effective Lagrangian ([^. 

The decays to one lepton and two quarks can be written: 

dV^N^l+ud) (TUnY 0--^-yi+yu) Ua , ^ rc n , ^2 

& = ^ (1-x + 9,)3 + + 

+ 120 : 2(2 - x)(l - a: + yi)y/yiyu + 0 : 3 ( 20 ;^ - a;^(5 + 5yi + yu) - %;(! + yi + 2yu) 

+ x{3 + lOyi + 3yf + 3yu + 3yiyu))] + 24Q;4a;(l -x + yify/my^} (10) 

with 2y/^ < X < 1 + yi — yu, yi = j ; W = the coefficients oi ,..,4 take the 


expressions: 


^S2,iu ^S2,iu^S3,iu) ^iu,i 


«i - + 

f yw{l -x + yi-yw) , . 

O 2 — I Clsi/., NO I ^S3,i„OCVo,iu I 

V [l-x + yi- ywY + ywyvw ) 

yl^il-x + yi -yw) 




{l-x + yi- ywY + ywyv^ 


0:4 Ois2,in^V,iu^in,il 

»i‘h!/H-=(tt)'. »r„ = (S^) 


dT^N^^dd) i'mN\^x^ ^{1 - x){l - X - Aya) ^ 2 ^ ^ , o , o 

& =1^1—It --K, J3 + x(-5 + 2 . + 2 ,,)) 


1287r3 V A / 4 {1 - x^ 

® i^s2,id ~ ^S2,id^s3,id) (1 ~ a;)(l — a: — 2yY)'\ 


( 11 ) 


with 0 < a: < 1 - Ayd, yd = 

rriN r _ 

1287r3 V A y' 


dx 1287r3 V A y V (1 - a;) ^ yu)o^„,^l 


( 12 ) 

















with 0 < a: < 1 - 4?/„. 


And the purely leptonic decay: 


leptons) 

dx 



with 2^1 <x <l + yi-yv, yi 
take the expressions: 





a2 — 12asQ^^,aw,ii 


V-^ I i/t yhh ^ 

[l-x + yi- ywY + ywWw 


{l-x + yi- ywY + ywVvw 
{l-x + yi-yw 


P{x) = 2x^ - x^(5 + 5yi + yi>) - %;(! + yi + 2yi') + a;(3 + IQyi + 3|/f + ?,yi' + 2,yiyi') 
R{x) = {2-x){l-x + yi)y/y^. 


In the last expressions x = ‘^P^epton/'^N- 


Finally, in the considered mass range, the one-loop operators in the Lagrangian (|^, 
induce the decay of N to neutrino and photon: 



(14) 


This decay mode leads to an interesting phenomenology, as will be shown in the following 
sections. 

C. Bounds on the couplings 

Existent bounds on right-handed heavy Majorana neutrinos (often called “sterile”, as 
they are SU{2) singlets) are generally imposed on the parameters representing the mixing 
between them and the light left-handed ordinary neutrinos (“active”). Very recent reviews 
[ini[nii2E] summarize in general phenomenological approaches the existing experimental 
bounds, considering low scale minimal seesaw models, parameterized by a single heavy 
neutrino mass scale and a light-heavy mixing Uin, where I indicates the lepton flavor. 
The mentioned mixings are constrained experimentally by neutrinoless double beta decay, 
electroweak precision tests, low energy observables as rare lepton number violating (LNV) 
decays of mesons, peak searches in meson decays and beam dump experiments, as well as 
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direct collider searches involving Z decays. Also, previous analysis [I8l[27j refer in general 
to similar heavy neutrino-standard boson interaction structures, e.g.: 


Lw = -^h^Ui^PLNW^ + h.c. 
v2 


(16) 


Pz — 


2cw 


i^LYUmPiNZf, + h.c. 


(16) 


The effects of this modihcation on the weak currents are studied, as they lead to correspond¬ 
ing variations in the weak bosons decay rates and W and Z mediated processes involved in 
the existing experimental tests, specially in colliders [161 EHHSl] . 

In the effective Lagrangian framework we are studying, the heavy Majorana neutrino 
couples to the three fermion family flavors with couplings dependent on the new ultraviolet 
physics scale A and the constants where i labels the families and O the operators. 

The operators presented in ([^ lead to a term in the effective Lagrangian (|^ that can 


be compared to the interaction in (|15|), and a relation between the coupling and the 

(i) 2 

mixing Uin was derived in [19]: Ul^N — "^2 , while no operators lead to a term that can be 


directly related -with the same Lorentz-Dirac structure- to the interaction in (16) (nor at 
tree or one-loop level). Some terms in the Lagrangian (j^ contribute to the ZNv coupling. 


but as they are generated at one-loop level in the ultraviolet underlying theory, they are 
suppressed by a l/167r^ factor. 

In consequence, we take a conservative approach. In order to keep the analysis as simple 
as possible, but with the aim to put reliable bounds on our effective couplings, in this work 
we relate the mixing angle between light and heavy neutrinos {UeN, Ufj,N, UrN) with the 
couplings as f/ ~ where v corresponds to the vacuum expectation value: v = 250 GeV. 
As we will explain shortly, we consider two situations in which the different bounds applies 
to the couplings. 

Some of the considered operators contribute directly to the neutrinoless double beta decay 
(Oz//3/3-decay) and thus the corresponding coupling constants, involving the hrst fermion 
family i = 1, are restricted by strong bounds. We explicitly calculated the implications for 
the effective couplings in our Lagrangian. 

In a general way, the following effective interaction Hamiltonian can be considered: 


l-i = Geff uTd eTN + h.c. 


( 17 ) 







where T represents a general Lorentz-Dirac structure. Following the development presented 
in |33] and using the most stringent limit on the lifetime for Oz//?/3-decay > 2.1 x 10^® 

years obtained by the Gerda Collaboration [SB] we have obtained the following bounds on 


the Geff 


Geff < 7.8 X 10“® 


mjv 

lOOGeC 


1/2 


GeV 


-2 


(18) 


The lowest order contribution to Oz/yg^j-decay from the considered effective operators comes 
from those containing the W held and the 4-fermion operators with quarks u, d, the lepton 
e and the Majorana neutrino N. These operators contribute to the effective Hamiltonian 
(17), with Geff = Thus we can translate the limit coming from Geff on which, for 
A = 1 TeV, is 


^ X 10“^ 


rriN 
100 GeV) 


\ 1/2 


(19) 



(a) (b) 


FIG. 2: (a)The Branching ratios for the Majorana neutrino decay with coupling constants 
in the set A (solid lines) and set B (dashed lines). The unlabeled curves correspond to the 
decay N —)■ i/uu and N —)■ udd. (b) Total width for the two combinations of coupling 
constants dehned in the text (set A and set B) and A = 1 TeV. 

As will be explained in Sec |TT| the relevant Majorana neutrino mass range for considering 
this heavy neutral particle as a solution to the MiniBooNE anomaly is 400MeV < m^r < 
600MeV [23]. The experimental bounds for this mass values are exhaustively discussed 
in mi and references therein. Taking into account that the MiniBooNE experiment deals 
with muon-type neutrinos, we now discuss the bounds on the mixings, which are not 
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constrained by Oz/^/ 3 -decay, and are most restrictive than the existing ones for the third 
fermion family. 

As can be seen in na, the existing bonnds for U^jsf for 500 MeV come from beam 

dnmp experiments as NnTeV [37], CHARM II [3H| and BEBC |39|, rare lepton nnmber 
violating (LNV) meson decays at LHCb [iQ] and from colliders as those from DELPHI |3l] . 
In the case of the heavy Majorana nentrino with effective interactions we are considering, the 


clear dominance of the nentrino pins photon channel fonnd in (14) makes the beam dnmp 


and rare LNV experiments bonnds inapplicable, as this decay mode to invisible particles is 
not considered in those analysis, and can considerably alter the nnmber of events fonnd for 
N decays inside the detectors mi EE], 

In the light of this discnssion, we consider the bonnds from DELPHI following 
the treatment made in [18] ■ Iii our case, for only one heavy Majorana nentrino we have: 

= Uip^Ui/N and the allowed valnes for the mixings are of order: 


5 X 10 


( 20 ) 


For the Lepton-Flavor-Violating processes e.g. fi —>■ ey, fi —)■ eee and r —)■ eee, which 
are indnced by the qnantnm effect of the heavy nentrinos, we have very weak bonnds for 

rriM < mw [SI 1^1^ - 


Thus, the bound in (20) can be translated to the constants a, and we have for A = 1 
TeV 


^ bound ^00 

Ocoii < 2.3 


( 21 ) 


For completion we have explicitly calculated the bounds that can be inferred from the 
single Z —)• uN and pair Z ^ N N “excited” neutrino production searches at LEP [33]. The 
hrst process can be generated by one-loop level effective operators (|^ giving the terms in 
the Lagrangian ([^. As the one loop level couplings are supressed by the factor I/IOtt^, the 
corresponding bound for the couplings is absorbed by the (Ibvr^)^ multpliying 


the bounds, so that the collider (21) value is still more stringent. It is important to mention 
that other effective operators (4-fermion operators in (j^) contribute to the vN and NN 


production at LEP, but at the Z peak they give less restrictive bounds than the ones in 


(20). For the decay Z ^ N N, we have a direct contribution from the tree level operator 


10 






Onn<i„ giving 




( 22 ) 


A conservative limit for any mAT mass is Sr(Z — )■ NN)Br‘^{N — )■ i>(9)'y) < 5x10“^ [33]. This 
result is model-independent and holds for the production of a pair of heavy neutral objects 
decaying into a photon and a light invisible particle. For the low mjv values considered in this 
work, we can take Br{N —)■ v{9)^) ~ 1 and the corresponding bound is < 3.0 x 10“®, 


more restricting than the bound in ( 20 ), but not taken into account, as the corresponding 


operator does not contribute to the N decay. 

In order to simplify the discussion, for the numerical evaluation we only consider the 
two following situations. In the set we call A the couplings associated to the operators that 
contribute to the 0z//?/3-decay (C>ve 0 , OduNe, Oq^nl, OlnqiI and Oq^lci) for the hsrst family 
are restricted to the corresponding bound and the other constants are restricted to 

the bound determined by colliders o^cou'^- called B all the couplings are 

restricted to the bound which is the most stringent. For the 1-loop generated 

operators we consider the coupling constant as l/( 167 r^) times the corresponding tree level 
coupling: = Q;*’'“/( 167 r^). Thus, for the operators Odn^ Onw and OpTv, which 

contribute to Of(3(3 we have 

1 


vd) n/d) n/d) ^ bound 


(23) 


for fermions of the hrst family. For the remaining operators we take 


bound bound 
^Coll ? 


(24) 


in the sets A and B respectively. 

Figure shows the branching ratio as a function of the Majorana neutrino mass rriN- The 
decay is calculated for different values of the constants a^. We show the branching ratios 
for both sets A and B. It can be seen that, for low masses, the dominant channel is the 
decay of N to photon and neutrino. Figure [2^ shows the total decay width dependence on 
the mass for both coupling sets considered. 

Taking the values of the couplings ^ be equal for every family i, and also for every 
tree level coupling and taking the one-loop generated couplings as 0 ;^“^°°^ = 


In Fig. ^ we show the results for the Majorana neutrino decay presented in Sec. IB 
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we derived an approximated expression for the ratio between the widths r(A^ —)■ in 


(14) and T{N ^ I'^ud) in (10): 



(25) 


This limiting value explains the behavior found in Figj^for low Majorana neutrino masses, 
showing the neutrino plus photon decay channel is clearly dominating. This is an interesting 
fact since we have a new source of photons, leading to a very rich phenomenology discussed 
in the next section. 

II. APPLICATION TO NEUTRINO-RELATED QUESTIONS 

Searches for heavy neutrinos often rely on their possibility to decay to detectable particles. 
The interpretation of the corresponding results for such searches requires a model for the 
decay of the heavy neutrino. Several explanations to different kind of problems seem related 
to weakly interacting neutral particles, like new neutrinos. In particular the MiniBooNE [21] 
anomaly or the observation of sub-horizontal air-showers by Cherenkov telescope SHALON 
[23] have possible explanations by long lived neutral particles like the one studied in this 
work. 

The MiniBooNE experiment was built to search for —)■ z/g conversion, in order to 

conhrm or refute the previous results of LNSD, which were inconsistent with global neutrino 
oscillation data [12]. The MiniBooNE anomaly consists in an unexplained excess of low 
energy electron-like events in charge-current quasi-elastic electron neutrino events over the 
expected standard neutrino interactions mm- 

This excess of electron-like events could be caused by the decay of a heavy neutrino. 
This solution was proposed by Gninenko m in a model with sterile neutrino mixed with 
the standard neutrinos by a matrix U. He hnds that N with 


400MeU < rriN < 600MeU 
10 “^ < |f/^7v|2 <410“^ 

< tn < 10“®s 


(26) 


could explain the anomaly, as the excess of electron-like events in the beam could be 
caused by the decay of a heavy neutrino with a radiative dominant decay mode —)■ z/y 


12 








where the hnal photon would be converted into an e’''e“ pair with a small opening angle, 
indistinguishable from an electron in the detector. This is called a converted photon. 

The Gninenko analysis is based on the assumption that the heavy neutrino radiative decay 
is dominant. The effect of the mentioned strong radiative decay is the flux attenuation by N 
decay and then the decrease of the signal events in the detector. The consequences are less 


restrictive bounds on [IZl ES], as we explained in secTC The proposal is then that 

the excess of events observed by MiniBooNE could originate from converted photons and 
not from electrons. The future experiment MicroBooNE will provide a test to this proposal, 
as it will be able to separate photons from electrons or positrons |43j . 

In the context of the effective interactions considered in this work, one has to check if the 
N —)■ uA is the dominant decay, by comparing the decay of N to pions, which is the correct 
hadronic hnal state for the low masses studied here. We have found that the corresponding 
decay is mainly given by 


■^N^l+n- 


Gl 


(z) 2 ' 

OlwV 


Stt \ 2A2 




(1 


m 


l \2 


m„ 


m 


N 


m 


N 


1 + 


mr 


m 


Af / J 


X 


1 + 


mf 


o \ 2 2 

mi \ , mf 

- ^-k- 


m 


N 


m 


N 


m 


(27) 


N 


In the mass range proposed m we hnd that the ratio of the branching ratios for the 
different decay channels is Br{N —)■ l^7i)/Br{N —)■ z/(p)A) ~ 8 x 10“® and Br{N —)■ 
leptons)/Br{N —>■ u{i?)A) ~ 4 x 10“® thus conhrming the dominance of the radiative decay 
N uA. 

The heavy neutrino N could be directly produced by the in neutrino-nucleon reactions 
by the effective operators Oqunl, OlnqcI and Oq^Ld, with the subsequent decay and photon 
conversion as we show in Fig.([^. 

The excess of events is related to the relative magnitude between the Standard Model 
Neutral Current (SM NC) —)■ process and the effective NC-like N production 

—>■ NAf being Af a nucleon. For the effective operator we have a 4-fermion contribution 
with intensity a/A^ and for the SM-NC j (4m^). Then the amplitude ratio is /C = an^/2A^ 
and the N production is weighed by the factor Kf = (an^/2A^)^ relative to the SM-NC 
scattering. The constant plays the role of the mixing matrix in the Gninenko [2l] 


work, and then the value found in (26) is consistent with the allowed value by the 


collider bound of (20) 
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FIG. 3: Production process of extra neutrino N by effective interactions and their 

subsequent decay. 


The constraint for the lifetime of the heavy neutrino in (26) must also be fulhlled in order 


to consider the N effective radiative decay as an alternative explanation for the MiniBooNE 
anomaly. In Fig.(|^ we show the lifetime tn as a function of mjq for the sets A and B and 
for A = 1 TeV. In the case of A > 1 TeV the allowed region is upwards the curves. Thus, 
we can see a region compatible with set A where tat < 10“® s as in the solution proposed 
by Gninenko. 

As was previously mentioned, this kind of neutral particle which decays dominantly to 
neutrino and photon could be the explanation for several sub-horizontal events detected by 
the Gherenkov telescope SHALON as it was recently proposed in [23]. In the cited work 
the authors propose that the solution could be a neutral and then penetrating long-lived 
massive particle able to cross 1000 km of rock and decay within the 7 km of air in front 
of the telescope. In Fig. we show the decay length as a function of the heavy neutrino 
mass for different energies and couplings in the sets A and B. We can see that there is a 
region of the parameter space which could possibly explain the SHALON observations with 

a Idecay ~ 1000 km. 


To conclude, a few words about the detectability of this particle in colliders like the LHG. 

Searches for neutral long-lived particles as the heavy neutrino proposed by [21| have been 
studied in the context of t~ rare decays jH], where the authors propose to search for events 
with two vertices, featuring the production and decay of the unstable neutrino N. The 
use of displaced vertices has also been proposed to search for sterile neutrinos at the LHG 
gsiiizi, for N decaying to leptons and quarks or purely leptonically. Early displaced vertices 
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FIG. 4: N lifetime as a function of its mass for the sets of coupling constants labeled A 
and B for A = 1 TeV. For A > 1 TeV the values allowed for the lifetime correspond to the 
region upwards the curves. The horizontal solid line corresponds to the limit value found 
in the Gninenko solution for the MiniBooNE anomaly (tat < 10“®s). 

searches are reviewed in |45] . 

As we have shown, for the N masses considered in this work the dominant decay is 
the radiative A —)■ z/y channel, which can be observed by the signature of an isolated 
electromagnetic cluster together with missing transverse energy: 

7 + (28) 

where the photon originates in a displaced vertex. 

New physics searches involving such hnal states have been performed at the LHG |50] , 
and it has been suggested that this signal could be enhanced with the combined use of missing 
transverse monentum plus photons and displaced vertices searching techniques [SI1E2]. The 
use of this technique will allow to probe parts of the parameter space which are inaccessible 
by other methods. The use of displaced vertices has the advantage that for decay lengths of 
the order of, very roughly L G (10“^ — 1) m, there is little standard model background. We 
hnd that decay lengths as the above mentioned for masses between 1 — 30 GeV are possible 
in this model as we show in Figj^for the sets A and B. 
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(a) N Decay Length 


(b) N Decay Length (higher masses) 


FIG. 5: Decay length for different nentrino energies as a function of the neutrino mass for 
the coupling constant sets A (solid lines) and B (dashed lines) for A = 1 TeV. The 
energies, E = 10” GeV, vary from left to right with increasing n (1-5). For A > 1 TeV the 
decay length corresponds to the top right region from the curves. 


III. SUMMARY AND CONCLUSIONS 

We have calculated the decay widths and branching ratios for a relatively light heavy 
Majorana neutrino (with rriTv < in an effective approach, considering its possible 

decays to fermions, quarks and photons, focusing on a relatively low neutrino mass range. 
We hnd that for masses below approximately 30 GeV the dominant channel is the neutrino 
plus photon mode: N —)■ vA. With this decay mode in mind, we explored the plausibility of 
considering it as an explanation for the MiniBooNE and SHALON anomalies. We checked 
that in the effective model the radiative decay is dominant respect to the lepton plus pion 
mode, and leads to values of the effective couplings a which are consistent with the mixing 
value found by Gninenko [21] and with collider bounds |3l|. Also, we show that 

the Majorana neutrino lifetime also hts the limits in [2l|. This kind of weakly interacting 
long-lived particle has also been proposed as an explanation for sub-horizontal events in 
the SHALON telescope [23], and we hnd that the N decay length is compatible with the 
proposed explanation for part of our parameter space. This kind of particle could also 
be searched for in the LHG, with the use of the displaced vertices technique, with little 
standard model background. 
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